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Abstract 

In this article we derive rigorously amplitude equations for stochastic PDEs with 
quadratic nonlinearities, under the assumption that the noise acts only on the sta- 
ble modes and for an appropriate scaling between the distance from bifurcation 
and the strength of the noise. We show that, due to the presence of two distinct 
timescales in our system, the noise (which acts only on the fast modes) gets trans- 
mitted to the slow modes and, as a result, the amplitude equation contains both 
additive and multiplicative noise. 

As an application we study the case of the one dimensional Burgers equation 
forced by additive noise in the orthogonal subspace to its dominant modes. The 
theory developed in the present article thus allows to explain theoretically some 
recent numerical observations from [Rob03|. 

1 Introduction 

Stochastic Partial differential equations (SPDEs) with quadratic nonlinearities arise 
in various applications in physics. As examples we mention the use of the stochas- 
tic Burgers equation in the study of closure models for hydrodynamic turbulence 
HCY951 and the use of the stochastic Kuramoto-Sivashinsky equation or similar 
models HCB95 1 ILCM961 IBGR021 iRML+OOl for the modelling of surface phenom- 
ena. Very often SPDEs have two widely separated characteristic timescales and 
it is desirable to obtain a simplified equation which governs the evolution of the 
dominant modes of the SPDE and captures the dynamics of the infinite dimen- 
sional stochastic system at the slow timescale. The purpose of this paper is to 
derive rigorously such an amplitude equation for a quite general class of SPDEs 
with quadratic nonlinearities and, furthermore, to obtain sharp error estimates. 
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Consider, as a working example of the class of SPDEs that we will consider in 
this paper, the following variation on the Burgers equation 

dtu = d^u + ud x u + (1 + 7)« + acj) (1.1) 

subject to external forcing a<p and to Dirichlet boundary conditions on [0, tt]. Since 
we are working very far from the inviscid regime, the solutions to this equation in 
the absence of forcing would decay quickly to 0, were it not for the extra linear 
instability (1 + j)u. The constant 1 appearing in this term is taken equal to the 
Poincare constant for the Dirichlet Laplacian on [0, it] and is designed to render the 
first mode sin(x) linearly neutral. The constant 7 therefore describes the linearised 
behaviour of that mode. The aim of this article is to study the behaviour of solutions 
to (11.11 ) for small 7 over (large) timescales of order 7 . 

It is well known HBlo051ICH93l (see also MSch94|[Sch0ni for general results on 
unbounded domains) that in the absence of forcing, the solution to (11.11 ) is of the 
type 

u(t, x) = yfya(rft) sin(x) + 0( 7 ), (1.2) 
where the amplitude a solves the deterministic Landau equation 

dta = a — j^a 3 . 

If the forcing <j) is taken to be white in time (actually, any stochastic process with 
sufficiently good mixing properties would also do), then, provided that a = 0(7), 
the solution to dl.ll ) is still of the type dl.2| ), but a now solves a stochastic version 
of the Landau equation: 

d t a = a - j^a 3 + a£(t) , 

where £ is white noise in time and the constant a is proportional one the one hand 
to the ratio a/7 and on the other hand to the size of the projection of cj> onto the 
'slow' subspace spanned by the mode sin(x) BBlo051 . In particular, one gets a = 
if the projection of <fi onto that subspace vanishes. 

This naturally raises the question of the behaviour of solutions to (11.11 ) when the 
external forcing acts only on the orthogonal complement of the 'slow' subspace. 
Roberts [Rob03] considered for example noise acting only on the second mode 
s'm(2x). Using formal expansions relying on centre manifold type arguments, he 
derived a reduced model describing the amplitude of the dominant mode. More- 
over he demonstrated numerically that additive noise is capable of stabilising the 
dominant mode, i.e. the noise eliminates a small linear instability. In turns out that, 
in order to have a non-trivial effect on the limiting amplitude equation, the strength 
of the noise should be chosen to scale like ^7, i.e. a = 0(^/7). We show that in 
this case, one has (after integrating against smooth test functions) 



u(t,x) = y/^Ofjt) sin(x) + C(7 5/8 ) 



(1.3) 
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To be more precise, we have additional noise terms of order y/j on higher modes 
that average out when integrated against test functions, i.e. they are small in some 
appropriate weak (averaged) sense. 

The amplitude a solves a stochastic differential equation of Stratonovich type 

da = (1 + 6{)a dt - ±a 3 dt + y/6 2 + S 3 a 2 o dB(t) . (1.4) 

Here, the constants 5{, i = 1,2,3 are proportional to a 2 /"/, a A /^ 2 , and cr 2 /7 
respectively, with proportionality constants depending on the exact nature of the 
noise. The Wiener process B can be constructed explicitly from the external forc- 
ing (j), but unless 62 = it is not given by a simple rescaling. 

In the particular case, where cj)(x, t) = sin(2x)£(i) with £ a white noise, one 
has 5i = —-§pj, 82 = and £3 = Note that 5\ is negative, so that if a 2 > 887, 
the solution to (11.41 ) converges to almost surely. This explains the stabilisation 
effect observed in IRob03l. 

In this article, we justify rigorously expressions of the form (11.31 ) for PDEs of 
the form (11.11 ) and we obtain formulas for the coefficients in the amplitude equation 
(11.41) . Unlike [B1605 ] we are interested in the situation where the noise does not act 
on the slow degrees of freedom directly but gets transmitted to them through the 
nonlinear interaction with the fast degrees of freedom. From a technical point of 
view, one of the main novelties of this article is that it provides explicit error bounds 
on the difference between the solution of the original SPDE and the solution of 
the approximating amplitude equation; this is a key requirement in tackling the 
infinite dimensional problem. Thus, our result is stronger in that respect than weak 
convergence type results in the spirit of e.g. ||EK86, Kur73J. Furthermore, we 
provide an explicit coupling between the two solutions, which is not trivial in the 
sense that, unlike in the case where the noise acts on the slow variables directly, 
the white noise driving the resulting amplitude equation is not a simple rescaling 
of the noise driving the original equation. 

Finite dimensional SDEs with quadratic nonlinearities and two characteristic, 
widely separated, timescales were analysed systematically by Majda, Timofeyev 
and Vanden Eijnden in a series of papers [MTE01 , MTVE99 ]. The SDEs that were 
studied by these authors can be thought of as finite dimensional approximations 
of stochastic PDEs with quadratic nonlinearities of the form (11.11 ) (in fact, the au- 
thors consider finite dimensional approximations of deterministic PDEs and they 
introduce stochastic effects by replacing the quadratic self-interaction terms of the 
unresolved variables by an appropriate stochastic process). In these papers, tech- 
niques from the theory of singular perturbation theory for Markov processes were 
used to derive stochastic amplitude equations with additive and/or multiplicative 
noise, which can be either stable or unstable. The results obtained by formal mul- 
tiscale asymptotics can be in principle justified rigorously using the theorem of 
Kurtz [Kur73], see also [EK86, Thm 3.1, Ch. 12]. However, since these results 
lack explicit error estimates, it is not clear a priori whether they can be applied to 
the infinite dimensional situation that we study in this paper. 
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The rest of the paper is organised as follows. In Section |2] we state the assump- 
tions that we make, and present our main result. Sections [3] to [5] are devoted to 
the proof of our main theorem. In Section [6] we apply our theory to the stochastic 
Burgers equation. 

2 Notations, Assumptions and Main Result 

The main object of study of the present article is the following SPDE written in the 
form (cf. IPZ921 ): 

du = {-Lu + B(u, u) + ve 2 u) dt + eQ dW(t) . (2. 1) 

Throughout this paper, we make the following assumptions. 

Assumption 2.1. L is a nonnegative definite self-adjoint operator with compact 
resolvent in some real Hilbert space TL. 

Let || • || be the norm and (•, •) be the inner product in TL. We denote by {ek}^ =1 
and {\k}kLi an orthonormal basis of eigenfunctions and the corresponding (or- 
dered) eigenvalues. We will furthermore assume that 

Assumption 2.2. The kernel of L is one dimensional, i.e. Ai = and A2 > 0. 

We will use the notation M := spanjei}, and P c for the orthogonal projection 
P c : TL — ► M. Furthermore, P s := I — P c . Before stating our assumptions on the 
nonlinearity B, we introduce the following interpolation spaces. For a > 0, we 
will denote by 7i a the domain of L a l 2 endowed with the scalar product (u, v ) a = 
(u, (1 + L) a v) and the corresponding norm || • || Q . Furthermore, we identify H~ a 
with the dual of TL a with respect to the inner product in TL. With this notation at 
hand we can state our next assumption. 

Assumption 2.3. There exists a G (0, 2) and j3 G (a— 2, a] such that B(u, v) : TL a x 
TL a — > TL^ is a bounded symmetric bilinear map with 

P c B(e k ,e k ) = 0, (2.2) 

for every k > 1. 

Finally, we assume that the Wiener process driving equation (12.11 ) satisfies: 

Assumption 2.4. W is a cylindrical Wiener process on TL. The covariance opera- 
tor Q is symmetric, bounded, commutes with L, and satisfies 

(e 1> Qei) = 0. (2.3) 

Furthermore, Q 2 L a ~ l is trace class in TL, where the value of a is the same as in 
Assumption I2.il 
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Remark 2.5. The scaling in e in equations f l2.il ) and \2.4\ below is dictated by the 
symmetry assumptions \2.2\ and \23\ . If either of these assumptions were to fail, 
the scaling considered in this article would not yield a meaningful limit. 

Remark 2.6. We could easily allow for a deterministic forcing term ef acting on 
the fast modes. We omit this for simplicity of presentation. 

We are interested in studying the behaviour of small solutions to (12.11 ) on 
timescales of order s~ 2 . To this end, we define v through ev(e 2 t) = u(t), so 
that v is the solution to 

dv = {-£- 2 Lv + vv + e~ l B{v, v)) dt + e~ l Q dW(t) . (2.4) 

Note that we made an abuse of notation in that the Wiener process W appearing in 
(12.41 ) is actually a rescaled version of the one appearing in (12. II ). but it has the same 
distribution. Now we are ready to state the main result of this article. 

Theorem 2.7. Let L, B, and Q satisfy Assumptions 12.7142.41 Fix a terminal time 
T > 0, a number R, as well as constants p > and k > 0. Then there exists 
C > such that, for every < e < 1 and every solution v of (12. 41) with initial 
condition v$ G 7i a and ||«o||a — R> there exists a stopping time r and a Wiener 
process B such that 

E sup \\P c v(t) - a{t)ei\\ p a < Ce vl ^ R , P(r < T) < Ce p . 

te[0,r] 

Here, a(t) is the solution to the stochastic amplitude equation 

da(t) = (ua{t) - fja{tf) dt + \J a h + a a a(t) 2 dB , a(0) = (v Q , e x ) , (2.5) 

where the coefficients u, fj, o~ a and aj, are given by equations ([4. 71 ), ( 14. #1 ), f |5.2| ), 
respectively. 

Furthermore, the fast Ornstein-Uhlenbeck process z(t) solving 
dz = -e~ 2 Lzdt + dW(t), z(0) = P s v , 
satisfies Esup te[0 r ] \\P s v(t) — z(t)\\a < Ce p ~ K . 

Remark 2.8. In a weak norm in time (for example H^ 1 ) one can show that z(t) is 
well approximated by white in time and colored in space noise of order e. Formally, 
we can write 

Z (t) = e~ tL£ ~ 2 P s v + e- 1 / e- (t - s)Le ~ 2 QdW(s) 

Jo 

» eL~ l Qd t W . 

Of course, for small transient timescales of order 0(e 2 ) the initial value P s vq of 
z(t) has a contribution of order 0(1). Thus estimates of the error uniformly in time 
are out of reach. 
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Remark 2.9. A immediate corollary of our result is that, under the assumptions of 
Theorem \2.7\ we can write 

E sup \\u(t) - ea(e 2 t)ei - eR(t)\\ p Q < Ce^~ K , 

te[0,re- 2 ] 

where u(t) is the solution to (12.11 ) with u(0) = 0(e), a(t) is the solution to the 
amplitude equation (12.51 ) with ea(0) = (u(0), e\) and R(t) = z(e 2 t) is the solution 
to 

dR = -LR + QdW, eR(0) = P s u(0). 

The noise that appears in the equation for R is a rescaled version of the noise that 
appears in the equation for z. 

Let us discuss briefly the main steps in the proof of this result. We first decom- 
pose the solution of (12.11 ) into a slow and a fast part: 

v(t) = P c v(t) + P s v(t) =: x(t) + y(t) , (2.6) 

to obtain a system of SDEs for (x, y), equation (I3.ll ). Our next step is to apply Ito's 
formula to suitably chosen functions of x and y in order to eliminate the 0(1/ e) 
terms from (13 -lb - We furthermore show that we can replace the fast process y by 
an appropriate Ornstein-Uhlenbeck process z. In this way, we obtain an SDE for x 
that involves only x and the (infinite-dimensional) Ornstein-Uhlenbeck process z. 
This is done in Section [3l see Proposition 13.91 

A general averaging result (with error estimates) for deterministic integrals that 
involve monomials of the infinite dimensional OU process z, see Corollary 14.51 
enables us to eliminate or simplify various terms in the equation for (x, z) and to 
reduce the evolution of x to the integral equation 

x(t) = x (0) + v f x(s) ds-fj I (x(s)) 3 ds + M(t) + R(t) , (2.7) 
Jo Jo 

where R(t) = 0(e 1 l 2 ~ K ) (for arbitrary k > 0) and M(t) is a martingale whose 
quadratic variation has an explicit expression in terms of (x, z). (We made an 
abuse of notation here and wrote x 3 for what should really be (x, ei) 3 ei.) This is 
done in Section |4] 

The final step in the reduction procedure is to show that the martingale M(t) 
can be approximated (pathwise) by the stochastic integral 

M(t)= [ s/a h + (x a a\s)dB(s), 
Jo 

where B(t) is a suitable one dimensional Brownian motion and a is the solution 
to the amplitude equation (12.5I ). This is done in Sections [5. II and [BT2l We remark 
that, whereas the derivation of equation (12.71 ) is independent of the dimensionality 
of x(t), the third part of the proof is valid only in the case where the kernel of L 
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is one dimensional. This is the price we have to pay in order to obtain rigorous 
explicit error estimates on the validity of the amplitude equation. 

Let us comment briefly on the case dim(ker(L)) = k > 1 but finite. The 
technique employed in the proof of Theorem 12.7 1 would still apply to this problem 
to give weak convergence of the projection of the solution of (12.41 ) to Ylj=i a j(t) e j> 
where a(t) would satisfy a vector valued amplitude equation. It does not seem 
possible, however, to obtain pathwise convergence using our approach, since the 
time change employed in the proof of Lemma 15.11 works only in one dimension. 
Neither does it seem straightforward to modify the present proof in such a way 
that one can obtain explicit error estimates without using Lemma IBTTI In the case 
where the amplitude a(t) is a Brownian motion on R. k , one could consider one 
dimensional projections, as was done in [HP04]. It is not clear however how to 
adapt the argument used in that paper to the case where the amplitude a(t) is the 
solution of a general SDE. Furthermore, the error estimate obtained in [HP04] 
scales like e c l k for some appropriate small constant c. This is cleary not optimal. 

3 The Reduction to Finite Dimensions 

Let us fix a terminal time T and constants k > and p > 0. Note that these 
constants are not necessarily the same as the ones appearing in the statement of 
Theorem 15 .21 above, but can get 'worse' in the course of the proof. 
Note first that one has 

Lemma 3.1. Under assumptions 12.71 \2~3\ and \2Z4\ equation (12. 4\) has a unique local 
(mild) solution u in TC a for every x G 7i a , i.e. u has continuous paths in 7i a . 

Proof. This follows from an application of Picard's iteration scheme for the mild 
solution, see for example [PZ92]. One can check that the assumption [3 < a — 2, 
together with the continuity assumption on B( ■ , • ), imply that the solution map 
has the required contraction properties for sufficiently small time. The fact that the 
stochastic convolution takes values in TL a is a consequence of Assumption l2.4l □ 

Remark 3.2. Note that we do not rely on a dissipativity assumption of the under- 
lying SPDE H2.4\) . Thus we can only establish the existence of local solutions. The 
existence of solutions on a sufficiently long timescale will be shown later to follow 
from the dissipativity of the approximating equations. 

Substituting the decomposition (12.61) into (12.41) . we obtain the following system 
of equations 

dx = vxdt + 2e~ 1 P c B(x, y) dt + e^PcBiy, y) dt (IXTk ) 
d y = (v- s~ 2 L)y dt + e~ l P s B{x + y,x + y)dt + e~ l QdW{t) . (l3Tb ) 

Since Lemma 13.11 does not rule out the possibility of a finite time blow up in TL a 
for the quite general system (13.11) . we introduce the stopping time 

r* = T A inf{t > | \\v{t)\\ a > e~ K } . (3.2) 
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Note that Lemma 13.11 ensures that, for a fixed initial condition vq and for e suffi- 
ciently small, one has r* > almost surely. 
Let us fix now some notation. 

Definition 3.3. For a real-valued family of processes {X e (t)}t>o we sa y X £ = 
0(f £ ), if for every p > 1 there exists a constant C p such that 

E( sup \X £ (t)\ P ) < C P f! . (3.3) 
v te[0,r*) ' 

We say that X £ = 0(f £ ) (uniformly) on [0,T], if we can replace the stopping 
time t* by the constant time T in A3. 3\) . If X £ is a random variable independent 
of time, we use the same notation without supremum in time, i.e. X £ = 0(f £ ) if 
E\X £ \p < C p fl 

We use the notation X £ = 0(f~) if X £ = 0(f £ e~ K )for every k > 0. 

3.1 Approximation of the stable part by an Ornstein-Uhlenbeck process 

In this subsection we show that the 'fast' process y(t) is actually close to an 
Ornstein-Uhlenbeck process, at least up to time r*. We have the following result. 

Lemma 3.4. Let z(t) be the Af 1 " -valued process solving the SDE 

dz(t) = -e~ 2 Lz dt + e^QdWit), z(0) = y(0) . (3.4) 

Then one has \\y(-) — z(-)\\ a = C(e 1_ ). 

Proof. It follows from the mild formulation of (12.41) that 



y(t) = z(t) + - f e -(*-^" 2 N(x(s), y(s)) ds , (3.5) 
£ Jo 

where we have used the notation N(x, y) = evy + P s B(x + y, x + y). From the 
properties of L we deduce that there exist positive constants C, c such that 

e F s \\ H ^ na <^ Ce _ ct fort > 1 _ (3-6) 
Since on the other hand Assumption [23] implies that 



\\N(x,y)\\ p <C(l + \\x + y\\ a ) 2 , 

the claim follows from the definition of r* and the fact that the right hand side of 
(13.61 ) is integrable for j3 > a — 2. □ 

The above approximation result enables us to obtain estimates on the statis- 
tics of the stopping time r*. For this we will need an estimate on the Ornstein- 
Uhlenbeck process (13.41) . 
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Lemma 3.5. Suppose that Assumption \L4\ holds. Then there is a version of z which 
is almost surely H a -valued with continuous sample paths. Furthermore, for every 
kq > and every p > 0, there exists a constant C such that 

E( sup \\z(t)\\ p a ) < Ce- K0 . (3.7) 
te[0,T] 

Proof. It follows for example from the proof of [PZ92, Theorem 5.9]. □ 

An immediate corollary of Lemmas |3.4| and 1331 is that the process y(t) is 'al- 
most bounded' in e. More precisely: 

Corollary 3.6. Assume that the conditions of Lemma \3.4\ and Lemma \3.5\ hold. 
Then, for every kq > and every p > 0, there exists a constant C such that 

E( sup \\y(t)Z) < Ce^°. (3.8) 
te[0,r*] 

Proof. Follows from Lemma [34] equation (13 -7b . and the triangle inequality. □ 

Note that the value of kq appearing in the statement above can be chosen in- 
dependently of the value k appearing in the definition of r*. Thus, with high 
probability, the event r* < T is caused by x(t) getting too large. To be more 
precise: 

Corollary 3.7. Under Assumptions 12.71 and \2.3\ for every p > there exists a 
constant C such that 

P(r* < T) < F (\x(t*)\ > e~ K ) + Ce p for e G (0, 1). 

Proof. Follows from Corollary 13 .61 and the Chebvshev inequality. □ 

3.2 Elimination of the 0(± ) terms 

Let us first introduce some notation. Given a Hilbert space TL we denote by TL (g) s TL 
its symmetric tensor product. Similarly, we use the notation v\ ® S V2 = |(i>i<8>i>2+ 
i>2 ® v i) for the symmetric tensor product of two elements and (A <g) s B)(x (8) y) = 
^ (Ax <S> By + By <S> Ax) for the symmetric tensor product of two linear operators. 

Let us recall that the scalar product in the tensor product space Ti a Hp is 
given by (u\ (g> v\,u 2 <& v 2 ) a ,p '■= (ui, U2)a{vi, v 2 )/3- With a slight abuse of 
notation, we write (•, -} Q := (•, -} a ,a- Furthermore, we extend the bilinear form B 
to the tensor product space by B(u (8) v) = B(u, v). 

With this notation, one can check tha{j]: 

Lemma 3.8. The operator (I(8) S L) _1 is bounded from TL^^sTL^ to TL 1+ \ ® s 7i 1+ \ 
for any 76K. 

'Since L has a zero eigenvalue, one should interpret (/ ® s L)^ 1 and L" 1 as pseudo-inverses, 
where, for instance, L _1 = on the kernel N(L). We will only apply these two operators to 
elements in A/" x so that this is of no consequence. 
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Proof. It suffices to note that I ® s L is diagonal with eigenvalues (Xj + Afc)/2 in 
the basis ej ® s eu- Note that N{I S L) = span(ei ® s e{). □ 

Now we are ready to present the main result of this subsection. 

Proposition 3.9. Let x and z be as above and let Assumptions \2~l\\2.4\ hold. Then, 
there exists a process R = C(e 1_ ) such that, for every stopping time t with t < t* 
almost surely, one has 

x(t) = x(0) + v f xds + A [ P c B(P c B(x, z), L~ x z) ds 
Jo Jo 

+ 2 1 P c B(x,L- 1 P s B(x + z,x + z))ds + 2 / P c B(P c B(z, z), lT x z) ds 
Jo Jo 

+ f P c B(I® s Lr 1 (z® s QdW(t)) + 2 [ P c B(x,L^QdW(s)) 
Jo Jo 

+ / P C B(I ® s L)~ l {z ® s P s B(x + z, x + z)) ds + R(t) . (3.9) 
Jo 

An immediate corollary is 

Corollary 3.10. Under the assumptions of Proposition \3.9\ define the process 
xr(V) by xji(t) = x(t) — R(t). Then for every p > and every a < 1/2, one 
has 

\xR(t) - x R (s)\ 
SU P ri r~ = 0(e u ) . 

0<s<t<T* \t — S\ 



Proof. It follows immediately from (13.91 ). using the definition of r*. The condition 
a < 1/2 arises from the two stochastic integrals in the right hand side of (I3.9I ). □ 

Proof of Proposition \3.9\ Applying Ito's formula to B(x, L^y), we get the fol- 
lowing identity in Hr\ 

dB(x, L~ x y) = 2uB(x, L~ x y) dt + 2e' 1 B(P c B(x, y), L~ x y) dt 
+ e^BiPcBiy, y), L~ l y) dt - e~ 2 B(x, y) dt 
+ e~ 1 B(x, L^PsBix + y,x + y)) dt + e~ 1 B{x, L~ l Q dW(t)) . 

Combining this with Lemma 13.41 and the continuity properties of B stated in As- 
sumption 12.31 it follows that, for every stopping time t with t < r* almost surely, 
one has 

2 I B(x,y)ds = 4e [ B(P c B(x,z),L~ 1 z)ds + 2e [ B(x, LT X Q dW(s)) 
Jo Jo Jo 

+ 2e B(x,L- 1 P s B(x + z,x + z))ds 
Jo 

+ 2e I B(P c B(z,z),L- 1 z)ds + R^t) , (3.10) 
Jo 
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where R^t) = 0(e 2 -). 

Applying Ito's formula to ~(y (£> y), we get the following identity in 7l a ~ 2 ® 

\d(y <S>y) = v(y ® y ) dt - e~ 2 (y ® s Ly) dt 

+ e~ l y ® s P s B(x + y,x + y)dt + e~ l y ® s Q dW(t) 

oo 

+ e- 2 ^Qei®Qeidt . (3.11) 
i=i 

Note however that all terms but the second one actually belong to TL^ 1 (g> TL a ~ l . 
Since B is bounded from TL a 1~L a into "HP and since (/ ® s L)^ 1 is bounded from 

n a-l jg n a-l tQ n a ^ w£ can aply p^j ^ jr j-1 tQ both s j des of fej^ 

Noting that P c B(ei (g) ej) = by Assumption 12.31 we get 

f P c B(y,y)ds = e [ P C B(I ® s Ly 1 ^ ® s P s B(x + z,x + z)) ds 
Jo Jo 

+ e [ P C B(I ® s L)-\z ® s QdW(t)) + Ra(t) , 
Jo 

where R2H) = 0{e 2 ~). Collecting both terms and inserting them into (13. lab con- 
cludes the proof. □ 

4 Averaging Over the Fast Ornstein-Uhlenbeck Process 

In this section, we simplify the equation for x further by showing that one can 
eliminate all terms in ( 13.91 ) that contain odd powers of z. Furthermore, concerning 
the terms that are quadratic in z, there exists a constant Q G TL a ® s TL a so that one 
can make the formal substitution z (g) z 1— ► Q. 

We start with a number of bounds on the integrals of products of Ornstein- 
Uhlenbeck processes. 

4.1 An averaging result with explicit error bounds 

Recall that QW can (at least on a formal level) be written as J2 < k'=2 Qk^kWkit) for 
some independent standard Wiener processes w^. For e > and k > 1, we define 
Zk{t) to be the stationary solution of 

dz k = -\k£~ 2 z k dt + q k £~ l dw k (t) . 

This is a Gaussian process with covariance 

E (z k (s)z k (t)) = . (4. 1) 

l\k 



Since z k (t) fluctuates very rapidly, one would expect from the law of large numbers 

that as e — > 0, one has z k — > 
precise by the following bounds: 



that as e — > 0, one has z k — ► weakly, but {z k ) 2 — ► weakly. This is made 
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Lemma 4.1. For every p > there exists a constant C p such that, for every t > 
s > and every k,£,m > 1, the bounds 



2 r ,2 r ,2 



E 



2Afc / / VAfcA^/ 
( f z k {r)z l {r)z m {r)dr) 2P < Cp [f^-)\t - sfe 2 ? 



- AfcA^A. 

/zo/<i. Here we denoted by 5 k i the Kronecker symbol. 

Proof. The first bound can be checked explicitly in the case p = 1 by using (14. II ). 
The case p > 1 follows immediately from the fact that f z k (r) dr is Gaussian. 

In order to obtain the other bounds, we recall first the fact that for an M 2p -valued 
Gaussian random variable X = (X\, . . . , X2 P ), we have 

p 

EX 1 -...-X n = II EX — X °* 

o-e£(2p) i=l 

where T>(2p) is the set of all permutations of {1, ... , 2p}. 

Turning to the second claim, consider first the case where k ^ I, so that z k and 
zi are independent. Thus 

2p 2p 2p 

E 11 h(ti)h(ti) = E f[ z k (ti)E f[ zi(ti) 

i=l i=l i=l 

^(fgy e n«p(->.-w, 

K ' crGS(2p)i=l 

where c = Ai is the smallest non-zero eigenvalue of L. The bound then follows by 
integrating over ti, . . . , *2p and using the fact that f* J s ' exp(— c\r — u\/e 2 ) dr du < 
Ce 2 (t-s). 

In the case where k = I, we have 

«=1 K Ac{l,...,2p} K i&A 

2Afc- 



E y-^t) E n E5 *^ s *^i-i)' 

Ac{l,...,2p} fc t6£ 2 (A)«=1 



where the sum runs over Y, 2 (A), the space of all permutations of numbers in A, 
where each number is allowed to appear twice. Now it is possible to check that all 
terms in the double sum where \A\ < 2p are cancelled by a term with a larger A, 
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where t T2i = t T2i _ 1 for some i. All remaining terms have |A| = 2p. It follows from 
(14.11) that there exists a constant C such that 



\t ai -ti +••• + !* 



o-GE(2p) 



The bound then follows immediately by integrating over ti, . . . ,t p . The last term 
can be bounded in a similar way. □ 

Let now Q £ H a ® s H a be given by 



00 ^2 



k=2 



2\u 



(4.2) 



The fact that it does indeed belong to H a ® s H a is a consequence of Assump- 
tion !2.4l Writing z(t) = Y^kLi %(*) e fc» we nave the following corollary of Lemma l4.ll 

Corollary 4.2. For every p > there exists a constant C p such that the bounds 



E 



z(r) dr 



E 



E 



z(r) (g) z(r) — Q) dr 
z(r) <8) z(r) <g> z{r)) dr 



2p 

a 


< C p (t - 


s) p e 


2p 




s) p e 


a 


< C p (t - 


2p 






a 


< C p (t - 


sfe 



hold for every t > s > 0. 

Proof. It is a straightforward consequence of the following fact. Let {vk} be a 
sequence of real- valued random variables such that there exists a sequence {7^} 
and, for every p > 1, a constant C p such that E|?jfc| p < C P 7fc- Then, for every 
p > 1 there exists a constant C' p such that 



E(TxtvlY<c'(yx 



/t=i 



fc=l 



*7i 



The result now follows immediately from Lemma 14.11 and from Assumption 12.41 
which states that the sequence q^X^ 1 is summable. □ 

Lemma 4.3. Let G £ be a famUy of processes in some Hilbert space K, such that, 
for every p > 1 and every k > there exists a constant C such that 



E 



G e (r) dr 



2p 



< C(t - sfe 2p , 



(4.3) 



Averaging Over the Fast Ornstein-Uhlenbeck Process 



14 



holds for every < s < t < 1. Then, for every p > and every k > 0, f/zere ex/sfs 
a constant C such that 



E( sup 



(n+l)5 



G £ (s)ds 



2p 



holds for every N > 0, every (5 € (0, (iV +1) ), a«<i every e > 0. 

Proof. It follows from ( 14.31 ) that, for every g > 0, there exists a constant C such 
that 

(n+l)<5 N 



P sup 



G £ (s)ds 



> K < CiV— , 

' - K q 



(4.4) 



holds for every K > 0. Note now that if a positive random variable X satisfies 
¥(X > x) < C/x q for every x > 0, then, for p < q, one has 



EX p = p / a^PCX > x) dx < 
Jo 



p / x^" 1 dx + Cp / xP- 9 " 1 



dx 



Combining this with (14.41 ) and choosing q sufficiently large yields the required 
bound. □ 

Proposition 4.4. Let K, he a Hilbert space, let f be a KL-valued random process 
with almost surely a-Holder continuous trajectories, let G £ be a family ofKL-valued 
processes satisfying (14. 3i . and let 



F £ (t):= [ (G £ (s), f(s)} ds . 
Jo 



Assume furthermore that, for every k > and every p > 0, there exists a constant 
C such that 

E sup \\G £ (t)\\ p < Ce~ K . (4.5) 
te[0,i] 

Then, for every 7 < 2a/(l + 2a), there exists a constant C depending only on p 
and 7 such that 

E||F £ ||P 1 _ a <C7(E||/||g i ) 1 /VP, 
where \\ ■ \\ c & denotes the a-Hblder norm for KL-valued functions on [0, 1]. 

Note that, if we can choose a < 1/2, but arbitrarily close, then we can choose 
7 < 1/2, but arbitrarily close, too. 

Proof. We focus only on the Holder part of the norm. The L°° part follows easily, 
as F £ (0) = 0, e.g. by taking s = in the following. 
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Choose 5 > to be fixed later. Moreover, for every pair s and t in [0, 1], set 
s = 5[5~ l s] and i = 5[5~ l t]. We furthermore define fg(t) = fit). One then has 



{G £ (r),f(r))dr < / (G e (r), (/(r) - f s (r))) dr + / (G e (r), f s (r)) dr 



<\\G E \\ L o 0(K) 6 & \t-s\\\f\\ c a+ / (G £ (r),f s (r))dr 



The second term in the right hand side can be bounded by 



\t-a\>6 



I (G £ (r),f s (r))dr + mm{\t - s\,25}\\f\\ cS \\G £ \\ L ~ (K) 

J s 



The first term of this expression is in turn bounded by 



2\t - s\ 



sup 

id- 1 



(n+l)(5 



G £ (r) dr 



nS 



C & 



Collecting all these expressions yields 



|F £ || cl _ a ^CWG.U^WfWc^ + d-'l sup 

V n<5" 1 



(n+l)6 



G £ (r) dr 



Choosing 5 = s 2 /^ 1+2a \ applying Lemma |431 and using ( 14.5b easily concludes the 
proof. □ 

We are actually going to use the following corollary of Proposition 



Corollary 4.5. Let z be as above and let a. be as in Assumptions \2.3\ and 12. 4\ Fix 

T > and let fi with i £ {1,2,3} be a-Holder continuous functions on [0,T] 
with values in {(J-i a )® 1 )*, respectively. Let F £ be given by 



F £ {t) ■= J Q {(fi(s))(z) + (f 2 (s))(z ®z-Q) + (f 3 (s))(z ®z®z))ds. 

Then, for every 7 < 2a/(l + 2a), there exists a constant C depending only on p 
and 7 such that 



E sup |i J, e (t)| p <C7e'n'(E(||/i|| ca + ||/ 2 || oa + ||/3|| <7fi 



,2ml/2 



where \\-\\c & denotes the a-Holder norm for ((7i a )® 1 )* -valued functions on [0, r*]. 

Proof. Note that z satisfies ( 14.51 ) with K, = TL a . This follows for example from 
the proof of [PZ92, Theorem 5.9]. The statement is then a consequence of Corol- 
lary 14.21 and of Proposition 14.41 □ 
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4.2 The reduction of the slow modes 

We now use the results of the previous subsection in order to show that most of the 
terms that appear on the right hand side of eqn. (13.91) are of order 0(e l / 2 ~). Note 
first that we can replace all occurrences of z in ( 13.91 ) by the stationary process z 
without changing the order of magnitude of the remaining term R. We are now 
ready to state the main result of this section. 

Proposition 4.6. Under Assumptions \2.1W2.4\ and with x and z defined as above, 
we obtain 

x(t) = x(0) + I P C B(I ® s L)- l {z{s) ® a Q dW{s)) 
Jo 

+ 2 P c B(x(s),L- 1 QdW(s)) 
Jo 

+ D j x(s)ds-fj (ei,x(s)fe 1 ds + R(t) , (4.6) 
Jo Jo 

where \\R(-)\\ = ©(e 1 / 2- ) and the constants v and fj are defined as 

D = u + 2(ei, B((I S L)-\B S ® s J) + (I ® L^B,) 

+ 2(B C L- 1 ))(ei ® Q)) , 
fj = -2(e 1 ,B(ei,L- 1 B s (ei,e 1 ))) . 

Here, we used the notation B s := P S B and B c := P C B. 

Proof. First we replace all instances of z by z on the right hand side of eqn. (13 -9b . 
which results in an error of order C(e 1_ ) which is absorbed into R. This is a 
straightforward but somewhat lengthy calculation which we do not reproduce here. 
We rely on 

z{t) - z{t) = e- tL£ '\P s v - z(0)) = e- tL£ ~ 2 0(e -) . 

Note that obviously \\z — z\\ ^ C(e 1_ ), as bounds uniformly in time are not avail- 
able due to transient effects on timescales smaller than 0(e 2 ). Nevertheless, we 
only bound the error in integrated form, which is sufficient for our application. Ac- 
tually, it is possible to check that the error terms which result from this substitution 
are of 0(e 2 ~). The only exception to this is the stochastic integral, where we ap- 
ply the Burkholder-Davies-Gundy inequality in order to get a remainder term of 
0{e 1 -). 

Once this substitution has been performed, the proposition follows from an 
application of Corollary !4.5l to the modified eqn. (13 -9b - The fact that, for every a > 
1/2, the various integrands are indeed d-Holder continuous functions with values 
in ((W a )®')* and Holder norm of order 0(e Q ~) follows from Corollary EM □ 



(4.7) 
(4.8) 
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5 Approximation of the Martingale Term 

This section deals with the final reduction step for the general system (13 -lb - We 
start by eliminating the stochastic integral of the type L z (g) QdW(s) from (I4.6I ). 
In fact, we show that we can replace the martingale part in eqn. (14.61 ) by a single 
stochastic integral of the type 



f V°a~+ 
JO 



abd 2 (s) dB(s) 



against a one-dimensional Wiener process B. Note that this section is superfluous 
in the particular case where the first stochastic integral in (14.61 ) vanishes. This is 
the case for example in the situation considered in [Rob03]. See Theorem 16. II in 
the next section. 

We emphasise that although all the previous steps are easily extended to higher 
dimensions, this step is valid only under the assumption that the kernel of L is 
one-dimensional. 

5.1 An abstract martingale approximation result 

We start with the following lemma; we will use it to approximate the martingale 
part of equation (14.61) by a stochastic integral against a one dimensional Brownian 
motion. 

Lemma 5.1. Let Mit) be a continuous Tt-martingale with quadratic variation f 
and let g be an arbitrary Tt-adapted increasing process with g(0) = 0. Then, 
there exists a filtration T% with Tt C Tt and a continuous Tt-martingale M(t) 
with quadratic variation g such that, for every 7 < 1/2 there exists a constant C 
with 

E sup \M(t) - M(t)\P < C(E 5 (T) 2 f) 1/4 (E sup \f(t) - g(t)\y 

t€[0,T] te[0,T] 

+ CE sup \f(t) - g(t)\P/ 2 . 
mo,T] 

Proof. Define the adapted increasing process h by 

h(t) = f(t) + miigis) - fis)) . 

s<t 

note that one has h < g almost surely. Furthermore, one has 

< hit) - his) < fit) - fis) 

for every t > s, so that one has < ^ < 1 almost surely. Define a martingale 
M(t) with quadratic variation h by the Ito integral 




!/(/! - / J^(s)dMis). 
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Define now an increasing sequence of random times Tt by 

T t = inf{s > | h(s) > g(t)} > t . 

Note that since h < g almost surely, the times T t are actually stopping times with 
respect to T, so that the time-changed process M (t) = M(T t ) is a martingale with 
quadratic variation g. Note that M (t) is a martingale with respect to the filtration 
T induced by the stopping times T t . Note also that T C T as a consequence of 
the fact that T > t almost surely. 

It remains to show that M satisfies the required bound. Let us start by defining 
the martingale A as the difference A = M — M. The quadratic variation (A) of 
A is then bounded by 



(a)w = i H^J dm 

ft / dh \ 
< J (l- ^r(s)) df(s) = f{t) - h(t) 

= sup(/0) - g(s)) . 

s<t 

Applying the Burkholder-Davies-Gundy inequalities [RY99J Cor. IV. (4. 2)] to this 
bound yields the existence of a universal constant C such that 

E sup |A(t)| p <CE sup \f(t) - g(t)\ p/2 . 
te[0,T] te[0,T] 

Before we turn to bounding the difference between M and M, we show that if 
F is an arbitrary positive random variable, B is a Brownian motion, 7 < 1/2, and 
q > p > 1, then there exists a constant C depending only on p, q and 7, such that 

E\\B\\ P ^ F < C{EF q ) p/2q , (5.1) 

where we defined 

\B(t) - B(s)\ 
\\ B h,F = sup — — — . 

0<s<t<F \ l — s \ ' 

One has indeed for every K > and every L > the bound 

P(||a|| 7l j > K) < F(F >L) + P(||5|| 7iL > K) . 

Applying Chebyshev's inequality and using the Brownian scaling together with 
the fact that the 7-Holder norm of a Brownian motion on [0, 1] has moments of all 
orders, this yields for every q > the existence of a constant C such that 

n\\B\\ 1>F >K)< inf ^- + " S£ ) < C- 



L>o\ L? K 2 i ) - Ki 
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The bound (15.11 ) then follows immediately from the fact that if a positive random 
variable X satisfies F(X > K) < (a/K) q for some a, some q and every K > 
then, for every p < q, there exists a constant C such that E|X| P < CdP . 

Note now that it follows from our construction that there exists a Brownian 
motion B such that M(t) = B(h(t)) and M(t) = B(g(t)). Noting that h < g and 
setting G = g(T), we have 

E sup \M(t)-M(t)\ p <E\\B\\ P sup \h(t) - g{t)\~< p 
te[0,T] ie[0,T] 

<E\\B\\ P G sup \f(t)-g(t)r 
te[o,T] 

and the result follows from (15.11 ) and Young's inequality. □ 
5.2 Application to the SPDE 

Before we state the next result, we introduce some notation. Let 7 G TL and 
T : H a -» H be defined by 

(y, 7) = 2(e 1 ,B(e 1 ,L~ 1 Qy)) , (y, Tz) = (e h B(I ® s L)~\z ® Qy)) . 



The facts that 7 € H and T is bounded follow from Lemma [X8] together with the 
fact that Assumption I2.4l implies in particular that Q is a bounded operator from TL 
to 7Y a_1 . Note that T is actually bounded as an operator from TL a ~ l to TL, but we 
will not need this fact. 

Theorem 5.2. Let Assumptions IX7T42. 3\ hold and let (x(t), y(t)) be the solution of 
(13.1b . Let z>, fj be given by (14.71) and (14.81) . respectively, and define 

o a = \\lf , °b = tr(TQr) , (5.2) 

where we identify Q from A4.2\l with the corresponding operato^from (TL a )* to 
TL a . 

Define finally X(t) = {x(t), ei). Then, there exists a Brownian motion B such 
that, if a is the solution to the SDE 



da(t) = ua(t) - fja 3 (t) + \fa h + a a a 2 (t) dB(t), a(0) = X(0), (5.3) 
then, for every p > and every k > 0, there exists a constant C such that 
E sup \X(t) - a(t)\ p < Ce p/i ~ K , 

for every e G (0, 1), where r* is defined in ( 13.21 ). 



2 An element u ® v of H a ® H a defines an operator from (H a )* to H a by (u ® v)(f) = u{f, v) 
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Proof. From Proposition 14. 61 we have that, with the notations introduced above, 

X(t) = X(0) + 9 [ X(s)ds- fj f X 3 (s)ds 
Jo Jo 

+ l (Tz(s),dW(s))) + [ X(s)( 7 ,dW(s)) + R 2 (t), 
Jo Jo 

where R 2 = 0(e x l 2 ~ K ). Denote by M(t) the martingale 

M(t)= [ (Tz(s),dW(s))} + [ X(s)( 7 ,dW(s)) . 
Jo Jo 

Its quadratic variation is given by 

f(t)= f \\ 1 X(s) + Tz(s)\\ 2 ds . (5.4) 
Jo 

It now follows from Corollary I4.5l that 

I/O - 9(-)\ = 0(e 1/2 ") , where g(t) = [ (a a X\s) + a h ) ds . (5.5) 

Jo 

Denote by M(t) the martingale with quadratic variation g(t) given by Lemma |5TTl 
and by x the solution to 

dx = vxdt- f]x 3 dt + dM(t) , x(0) = x(0) . 

It follows from Lemma O that M(t) - M(t) = 0(e 1/4 ~). Therefore, using a 
standard estimate stated below in Lemma 1531 

X (t) - x(t) = 0(e 1/4 ") . (5.6) 

The martingale representation theorem [RY99, Thm V.3.9] ensures that one can 
enlarge the original probability space in such a way that there exists a filtration Tt, 
and an JT t -Brownian motion Bit), such that both x(t) and x(t) are ^-adapted and 
such that 

dx{t) = ux(t) dt - rjx 3 (t) dt + V ' <*b + o a X\t) dB(t) . 

Note that in general the cr-algebra T% is strictly larger than the one generated by 
the Wiener process W up to time t. This is a consequence of the construction of 
Lemma IBTTl where one has to 'look into the future' in order to construct M. 
We finally define the process a as the solution to the SDE 



da(t) = Da(t) dt - fja 3 (t) dt + yfa b + a a a 2 (t) dB{t) . 
Denote p = a — x and G(x) = y&b + o a x 2 . Then, one has 
dp 2 (t) < 2vp 2 {t) dt + (G(a(t)) - G(X(t))) 2 dt + 2p(G(a(t)) - G(X(t))) dB(t) . 
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Using the fact that G is globally Lipschitz, this yields the existence of a constant 
C such that 

dp 2 (t) < Cp 2 (t) dt + C\X(t) - x(t)\ 2 dt + 2p(G(a(t)) - G(X(t))) dB(t) . 

It is now easy to verify, using Ito's formula, (15.61 ), and the Burkholder-Davies- 
Gundy inequality, that 

p = 0{£ l/i ~) and thus X(t) - a(t) = 0(£ l/i ~) , 

which is the required result. □ 

Let us finally state the a-priori estimate used in the previous proof. 

Lemma 5.3. Fix v E M and r] > 0. Let Mi(t) be martingales (not necessary with 
respect to the same filtration), and Xi(t), i = 1, 2 be solution of the following SDEs 

dxi(t) = vxi(t) dt - T)Xi(t) dt + dMi(t) , (5.7) 

with xi(0) — X2(0) = ©(e 1 / 4- ) and xi(0) = O(e ~). Suppose furthermore 
Mi(t) = 0(£°-) andM x (t) - M 2 (t) = OCe 1 / 4 "), then 

xi(t) - x 2 (t) = Oie 1 ^-) . 

Proof. This is a straightforward a priori estimate which relies on the stable cubic 
nonlinearity in (15.7I ). First, one easily sees from Ito's formula that Xi(t) = O(e ~). 
Then using the transformation Xi(t) = Xi(t) — Mi(t) to random ODEs for xi(t), 
we can write down an ODE for the difference x\{t) — X2(t), which we can bound 
pathwise by direct a priori estimates. We will omit the details. □ 

5.3 Main Result 

Let us finally put the results obtained in this and the previous two sections together 
to obtain our final result for the system of SDEs (13.11) . 

Theorem 5.4. Let Assumptions \2. 1W2.4\ be true. Let (x(t), y(t)) be a solution of 
dj.il) . Furthermore, let z(t) be the OU -process defined in (13. 41 ) and r* the stopping 
time from ( IJ.2D . Let finally a a , a^, v, and fj be defined in ( 15.21 ). ( 14. 71 ). and ( 14.81 ) re- 
spectively. Then there exists a Brownian motion B(t) such that, ifa(t) is a solution 
of 

da(t) = ua(t) - f]a 3 (t) + \fu b + a a a 2 (t) dB(t) , a(0) = (x(0), ei) , (5.8) 

then for all T > 0, R > 0, p > and k > there is a constant C such that for all 
e G (0, 1) and all ||x(0)|| Q < R and ||y(0)|| a < R we have that 

P(r* > T) > 1 - Ce p , 

E sup \\x(t)-a(t)ei\\ p a <Ce p/i - K , and E sup \\y(t)-z(t)\\ p a < Ce p ^ K . 
te[0,r*] te[0,r*] 
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Proof. The approximation of y(t) by z(t) is already verified in Corollary 13.71 and 
Lemma 13.41 The approximation of x(t) by a(t) follows from Theorem 15.21 The 
bound on the stopping time r* follows then easily from the fact that a(t) = 0(1) 
and z(t) = 0(1) uniformly on [0, T]. □ 

Remark 5.5. With the notation B k e, m = (B(ek, e£), e m ), formulas d5.2| ), H4. 71 ), and 
(14. 8\) for the coefficients in the amplitude equation can be written in the form 

T, - „±f 2B lnll , B kll B eek q] ^ 2B M1 B kU qj ,™ 

^ x 2 4^ X k Xi ^ X k + X t X k ' 
k=2 k k,e=2 K 1 k,e=2 K 

E2B k nBnk ir7xL\ 
r , (E9b) 
Ah 

k=2 fc 
00 4R2 J 00 o R 2 2 2 

A * (Afc + A - } Afc 



0",, 



If one chooses to expand the solution in a basis which is not normalised, i.e. one 
takes basis vectors = Ck&k, then the coefficients appearing in the right-hand 
side of the equation for the expansion transform according to 

5 r> CfcQ „ % 

-Dfcfm — -Dfcfefi , % — — • 

Cm Cfc 

/f is straightforward to see that v and a a are unchanged under this transformation, 
whereas fj is mapped to c\f\ and 05 is mapped to o~\,jc\ as expected. 

6 Application: The Stochastic Burgers Equation 

In this section we apply our results to a modified stochastic Burgers equation: 

du = (dl + l)u dt + ud x u dt + e 2 vu dt + eQ dW, (6. 1) 

on the interval [0, n], with Dirichlet boundary conditions. We take 

H = L 2 ([0,7r]) , e k {x) = ^-sm(kx) : 

B(u, v) = ^d x (uv) , L = -8 2 x -l, X k = k 2 - 1 . 

We also take W(t) to be a cylindrical Wiener process on TL and Q a bounded oper- 
ator with Qei = and Qe^ = qk&k f° r k > 2. It follows that 7i a = Hq([0, it]) is 
the standard fractional Sobolev space defined by the Dirichlet Laplacian on [0, 7r]. 
With this choice, using the notation from Remark 1531 we get 

Bum = ~ + t\h+l,m ~\k- £\S\k-e\, m ) , (6.2) 
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where S^t is the Kronecker delta symbol. 

It is possible to check that assumptions 12.11 12.21 and 12.31 are satisfied for any 
a > since, for smooth functions u, v, and w, one has for example 

r p7r 

(uv)'(x) w(x) dx = \ u(x)v (x)w'(x) dx < \\u\\ 
Jo 

< C\\u\\ \\v\\ \\w\\-hi , 

provided that one takes 7 < —3/2. (Values of a other than can be obtained in a 
similar way by using different Sobolev embeddings, see also [DPDT94].) Whether 
the trace-class assumption on Q 2 L a ~ l is satisfied or not depends of course in a 
crucial way on the coefficients {q^'kLi- 

The following result justifies the formal asymptotic calculations presented in 
HRob03ll . 

Theorem 6.1. Let u be a continuous Hq([0, n])-valued solution of d6.il) with ini- 
tial condition u(0) = 0(e), and assume that the driving noise W is given by 
a sin(2x)w(t) for a standard one-dimensional Wiener process w. Then there are 
Brownian motions B(t) and (5{t) (not necessarily adapted to the same filtration) 
such that if a is the solution of 

da = (v \ a dt a 3 dt H — \a\ o dB, ea(0) = — (n(0), sin(-)) L 2 

\ 88 / 12 6 tt 



and 



R(t) = - e - Lt P s u(0) + ( / e- 3(t - s) d[3(s) ) sin(2-), 



then for all K,p > there is a constant C such that 



E( sup \\u(t) - ea(e 2 f)sin(-) - eR{t)\\ p m ) < Ce^' 

H£[0,Te- 2 ] ' 



Proof. Note first that Assumption 12.41 is true for all a is this case, so that all the 
assumptions of Theorem 15 .41 are satisfied. Furthermore, we can use formulas ( 15.91 ) 
to obtain H, 

1 a 2 a 2 a 2 

77= , On = , <T/, = 0, V = V -\ . 

' 12 ' 36 72 88 

Note that the second term in the expression for v gives the Ito-Stratonovich correc- 
tion. 

However, the claim does not follow immediately, since we wish to get an error 
estimate of order e 3 / 2 instead of e 5 / 4 . Retracing the proof of Theorem 15 .41 we see 



3 Notice that a is the amplitude of the mode sin(a;) which is not normalized. This is in order to be 
consistent with earlier works on the stochastic Burgers equation. The modification of the formulas 
for the constants that appear in the amplitude equation in this situation is given by Remark 1531 
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that the claim follows if we can show that \f — g\ = 0(e ), where / and g are as 
in (15.41 ) and (15.51) . In our particular case, one has 7 = 0, so that 



fit) 



\Fz(sWds 



g(t) = a b t . 



The result now follows from Lemma [6721 below. 



□ 



Lemma 6.2. Let z be as in Corollarv \4.2\ Then, for every final time T, every p > 
and every k > there exists a constant C such that 



E sup 

te[0,T] 



(z(r) £3 z(r) — Q) dr 



2p 



< Ce 2p - K 



(6.3) 



Proof. We subdivide the interval [0, T] into N subintervals of length T/N, and we 
use the notation = kT/N. Using exactly the same argument as in the proof of 
Lemma |4~3l we see that, for every p > and every k > there exists a constant C 
such that 



E sup 

ke{0,...,N} 



(z(r) <8> z(r) — Q) dr 



2p 



< CN K e 2p . 



On the other hand, we know from Lemma [3751 that the W Q -norm of the integrand 
in (16.31 ) is of order O(s ~) uniformly in time. The claim then follows by taking 

N^e- 1 . □ 

Remark 6.3. In the case where only the second mode is forced by noise, one can 
actually take (3{t) = B(t), and (3{t) could be chosen to be a rescaled version of the 
Brownian motion that appears in equation (I6.ll ). 

The following theorem covers the case where W(t) is space-time white noise 
which is constrained to be antisymmetric around x = it. This corresponds to the 
case qj. = a for all k > 2. It is easy to check that Assumption 12.41 is satisfied for 
all a < \. 



We again use Remark [5751 with = \Jtt/2 in order to compute the coefficients. 
We obtain 



>1 



1 

12 ' 



a 

187T ' 



a b = Cb(J 



Cb 



9^2 ' 



2tt 2 (2k 2 + 2k + l)(k 2 - l)(k 2 + 2k) 

K=2 



and finally, 



v — v 



a o" 
367r 4tt 



E 

k=2 



k-1 k(k + l)J 2k 2 + 2k-l ' 



Note again that the first term in this expression is the Stratonovitch correction for 
the multiplicative noise term. This finally leads to 
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Theorem 6.4. Assume that a G [0, omd let ube a continuous Hq(\Q, ir])-valued 
solution of < \6. 1 I ) with initial condition u(0) = 0(e). Assume furthermore that the 
covariance of the noise satisfies = a for k > 2. Then there is a Brownian 
motion B(t) (not necessarily adapted to the filtration ofW(t)) such that ifa(t) is a 
solution of 

da(t) = va(t) - fja 3 (t) + \J a h + a a a 2 {t) dB(t), ea(0) = -(u(0), sin(-)) r 2 

7T 

where the constants are defined above, and 

1 r l 

R(t) = - e - tL P s u(0) + / e- {t - s)L QdW(s) , 
e Jo 

then for all K,p > there is a constant C such that 

E( sup \\u(t) - ea(e 2 t)sin(-) - eR(t)\\ p Ha ) < Ce^~ K 

He[0,Te- 2 ] ' 

for e sufficiently small. 
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